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Abstract This work is devoted to the construction and 
study of a mathematical model of hydropower unit, con¬ 
sisting of synchronous generator, hydraulic turbine, and 
speed governor. It is motivated by the accident hap¬ 
pened on the Sayano-Shushenskaya hydropower plant 
in 2009 year. Parameters of the Sayano-Shushenskaya 
hydropower plant were used for modeling the system. 
Oscillations in zones, which were not recommended for 
operation, were found. The obtained results are consis¬ 
tent with the full-scale test results carried out for hy¬ 
dropower units of the Sayano-Shushenskaya hydropower 
plant in 1988. 

Keywords Sayano-Shushenskaya hydropower plant, 
hydropower unit, synchronous generator, hydraulic 
turbine, speed governor, oscillations 


1 Introduction 

Nowadays one of the most important source of elec¬ 
tricity is the electric energy produced by hydroelectric 
facilities. According to [41], Norway gets 99% of its 
electric power from water, Brazil 84%, Austria 59%, 
Canada 58% and Russia 18%. Failures of hydropower 
plants cause loss of lives, major damage in the sur¬ 
rounding area, and serious economic consequences. In 
recent years, accidents at the hydropower plants have 
become frequent (see, e.g., Bieudron Hydroelectric Power 
Station (Switzerland, 2000), Taum Sauk Hydroelectric 
Power Station (Missouri, USA, 2005), Sayano-Shushens¬ 
kaya Dam (Russia, 2009), Itaipu Dam (Brazil, 2009), 
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Srisailam Dam (India, 2013), Dhauliganga hydro elec¬ 
tric station (India, 2013)). In order to prevent such ac¬ 
cidents, it is necessary to investigate their causes. 

This work is motivated by the accident happened 
on the Sayano-Shushenskaya hydropower plant in 2009 
year. According to the act of special commission of the 
Russian Federal Environmental, Industrial and Nuelear 
Supervision Serviee, immediately before the accident 
the power of the second hydropower unit was 475 MW 
at a head of 212 meters [36], i.e., it worked in the not 
recommended zone H (Fig. Q. Zone H of hydropower 
unit work is characterized by strong hydraulic turbine 
blows in flowing part and vibrations. For the normal 
operation it is recommended power range, correspond¬ 
ing to the zone HI, in which the efficiency of turbines 
has a maximum value. Also operation is allowed in the 
zone I, in which the dynamics is allowed, but the level 
of efficiency of the turbines are low. Operation in zone 
IV is not allowed. These work zones of hydropower unit 
of the Sayano-Shushenskaya hydropower plant were ob¬ 
tained by the full-scale test of hydropower units in the 
late 80s of 20th century and were published in the tech¬ 
nical report ^^Full-seale testing of turbines of Sayano- 
Shushenskaya hydropower plant with standard runneF 
No. 1008 d^. 

It is known [36] , that the accident was caused by vi¬ 
brations in the hydropower unit. Various possible rea¬ 
sons of occurrence of such vibrations were discussed 
in [361IT8] . However main attention in these works was 
given to the dynamical processes in the turbines and 
the full mathematical model turbine-speed governor- 
synchronous generator was not considered. 

In this paper an adequate mathematical model turbine- 
speed governor-synchronous generator is suggested. This 
model joins together differential equations of turbine 
(Pervozvanskiy’s model) [35], differential equations of 
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Fig. 1 Operational zones of hydropower unit of the Sayano- 
Shushenskaya hydropower plant 


synchronous generator (Park-Gorev model) [HIS], and 
differential equations of speed governor [33|. Parame¬ 
ters of the model are taken from m- As a result of using 
numerical methods the periodic solutions of differential 
equations of the model were found during the analy¬ 
sis of transient processes. These solutions correspond 
to hydropower unit vibrations observed for some oper¬ 
ation modes at the Sayano-Shushenskaya hydropower 
plant. 

Remark, that we considered also some other more 
simple models of hydropower unit (based on more sim¬ 
ple models of turbines (see [33]), synchronous gener¬ 
ators (the Tricomi equation [^1^ . the equations of 
synchronous generator with parallel connection in feed 
system [25]). However, for these models the above ef¬ 
fects have not been found. 


2 Components of hydropower unit 

Following m, let us consider the main constructive el¬ 
ements of hydropower unit. The hydropower unit in 
Fig-i consists of synchronous generator and radial- 
axial hydraulic turbine. In this work a nonregulated 
synchronous generator is considered. The rotor of gen¬ 
erator and the runner are connected together by a rigid 
shaft. 

The dam creates a difference in water level between 
the upper reservoir and lower reservoir. The flow of 
water is delivered from upper reservoir to turbine by 
penstock through the spiral casing. Water jets impact 
on the blades of the turbine producing torque applied 
to the rotating shaft. Since the turbine shaft is rigidly 
connected with the generator rotor, the rotor starts to 
rotate and to produce electricity, which is transferred 
to the grid. The water flow is controlled by means of 
guide vanes. 

For the safety of the power network, the frequency 
should remain almost constant. This is reached by keep¬ 
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Fig. 2 Scheme of the Sayano-Shushenskaya hydropower 
plant: 1 - upper reservoir, 2 - lower reservoir, 3 - dam, 4 
- penstock, 5 - blade wheel, 6 - guide vanes, 7 - rotor of 
synchronous generator, 8 - shaft, 9 - control gate 

ing the same speed of the synchronous generator. The 
rotational speed is controlled by the speed governor. 

The main structural elements of hydropower unit 
are presented in Fig. Introduce the following nota¬ 
tions: 5 is a signal of angular (rotational) speed devi¬ 
ation, /i is a position of guide vanes, Mt is a turbine 
torque, Mq is a generator torque, 1/ is a voltage re¬ 
quired by the power network. The power network rep¬ 
resents a set of energy suppliers and consumers. 



Fig. 3 The main structural elements of hydropower unit 

Thus, the hydropower unit represents a system con¬ 
sisting of hydraulic turbine, which produces a mechan¬ 
ical torque, generator, which converts mechanical en¬ 
ergy into electrical energy, and automatic speed gov¬ 
ernor, which regulates the rotation speed of hydraulic 
turbine. 

In order to develop a mathematical model of hy¬ 
dropower unit it is necessary to describe each structural 
element of hydropower unit, presented in Fig. Let us 
consider each elements in details. 


3 Mathematical model of synchronous 
generator 

The main constructive elements of synchronous gener¬ 
ator are stationary stator (Fig. a) and rotating rotor 
(Fig. a b). The windings are placed in the stator and 
rotor slots. Stator windings are arranged in such a way 
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that in the case of alternate current they generates a 
rotating magnetic field. 




a b 

Fig. 4 Synchronous generator, (a) - rotor: 1 - field winding, 
2 - damper winding; (b) - stator: 3 - stator winding 


Using Park’s transformation, the three-phase wind¬ 
ings of stator can be substituted by two equivalent 
short-circuited windings, and the rotor can be described 
as two equivalent short-circuited damper windings and 
one field winding. Let us choose the rotating coordinate 
system (d, g), connected with the rotor. 

q 



Fig. 5 The rotating coordinate system (d, q): axes (rd, rq) 
are rigidly connected with rotating rotor; axes (d, q) are 
rigidly connected with rotating magnetic field of stator 


Introduce the following notations: 

= Oq 


uj{t) = CJo + 


s{t) 


LOo 


where 0 is an angle between rotor axes and rotating 
magnetic field, generated by the stator windings [rad], 
Oq is an operating electrical angle [rad] (0 < Oq < 7r/2), 
0/^ is a deviation of operating angle, uj is an angular 
speed of rotor (shaft) [rad/s], ujq is a rated angular 
speed (synchronous speed) [rad/s], is a deviation 
of rated angular speed, s is a relative deviation of rated 


angular speed [p.u.]. The rated angular speed is deter¬ 
mined by the formula 

27r/ 

^0 =-, 

P 

where / = 50 is the network frequency [Hz = s“^], p 
is the number of pairs of poles. Since hydrogenerators 
of the Sayano-Shushenskaya hydropower plant have 21 
pairs of poles, then ujq = 27rb0l21 = 27r • 2.38 [rad/s] or 
ujQ = 27r3000/21 = 27r • 142.85 [rad/min]. Synchronous 
generators are usually operated with the nominal power 
factor cos^o = 0,8...0,9 m- For generators of the 
Sayano-Shushenskaya hydropower plant Oq = arccos(0, 9) 
Thus, a synchronous generator can be described by 
the Park-Gorev equations in per unit values [U pp.l32- 
149], [3 pp.5.15- 5.17]: 


= -^o(l + s)^q -ujorid-ujo Ud, 
iFg = CJo(l + s)^d - i^oriq- iAo Uq, 

= -^{Er — Eq)^ 


1 

Trd 


( 1 ) 


Erd — rr Erq-) 


Erq — rp Ej-di 

1 j^q 

where flux linkages are determined as follows 
Ed — Xd id T Eq fi- E^q^ 

Eq = Xq iq Erpd-) 

% = ^id + Eg + —Erg, 


^rd = ^id +Erg + ^ (Eg + Erg) , 
^rd ^rd 

'Erg = ^ig " Erd - ^Erd- 


( 2 ) 


vrq 


urq 


Here the following variables and coefflcients are rel¬ 
ative to the corresponding base values (voltage, current, 
flux linkage, impedance, inductance, power): r is a sta¬ 
tor resistance, id^ iq are currents in stator windings, Ud^ 
Uq are stator voltages, ^di Eq^ E^d^ E^q are flux link¬ 
ages, Er is a fleld voltage, Eq^E^d^E^q are electromo¬ 
tive forces, induced in the stator by the magnetic fleld 
of rotor winding currents for synchronous rotor speed, 
Xd^ Xq are synchronous inductances (reactances) along 
the axes d and g, Xrd^ Xrq are impedances of fleld 
winding, damper windings along the axes d and g, Xad^ 
Xaq are impedances of stator winding along the axes d 
and g. The following coefflcients are time constants: 
is a field-winding time constant with open stator and 
damper windings [s], T^d^T^q are damper winding time 
constants with open stator and fleld windings [s]. 
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The stator voltages Ud and Uq along the d- and q- 
axes in per unit values are determined according to the 
laws 


Ud = -U sin(6>o + Oa), Uq = U cos(6>o + Oa), 

where ^ is a voltage relative to the base voltage 

[p.u.], is a voltage in power network [V], Uh is the 
base voltage [V]. 

The motion of synchronous generator rotor about 
shaft is described by the torque equation in physical 
values pp.132-149], [3 pp.5.15-5.17]: 

Oa = 

J Ld = Mt - (^^d iq — • 

Here J is a moment of inertia of rotor (is a moment 
of inertia of hydropower unit) [kg • m^], ^di ^q are flux 
linkages in physical unit values [Wb = ^^ 2 ^ ]; 'iq 
are currents in stator windings in physical unit values 
[A], Mt is the rotation torque ]. Since the shaft of 

turbine is rigidly connected with the rotor of generator, 
a moment of inertia of rotor coincides with a moment of 
inertia of hydropower unit and an angular rotor speed 
coincides with an angular turbine speed. In our case 
the rotation torque Mt is a turbine torque, which is 
created by the pressure of water on the runner. 

Rewrite system of equations ^ in terms of base 
variables. For this purpose the second equation is di¬ 
vided twice by the base voltage: 


Uh = 

where iFb is the base flux linkage [Wb = ^b 

is the base impedance (valid also for resistances and 
reactances) [Ohms], Jb is the base current [A]. Then 


J 






Mt _^ f ^d ^q 

LdO^\yZ\yI\y C^Q^b \ ^b ^b 


^b 4 y ’ 


J cjq . Mt 

4b4 ^b4 


{^diq-^qid)- 


Thus, as a result of given transformation one obtains 
system of equations ^ in per unit values: 


Oa = 


(4) 


where Tj = is inertial constant of hydropower 

unit [s], 5'b = 4b4 is the base power [W = ]. 


4 Mathematical model of hydraulic turbine 


Hydraulic turbines can be classified by their type of 
construction, the most important ones being the Fran¬ 
cis, Felton and Kaplan or Propeller turbines [41]. They 
are distinguished by construction of runner and control 
methods of the rotation speed of the turbine. However, 
regardless of the turbine structure, the rotation equa¬ 
tion of turbine shaft in earth-fixed coordinate system 
has the following form nain] 

J uj = Mt — Mq^ 


where Mq is a resistance torque, which is created by 
the generator. 

The rotational torque is defined by the relation 


Mt = k 


Prit) 

L0{t) ’ 


where Pt is the turbine mechanical power [W = ], 

/c is a constant depending on the construction of turbine 
[p.u.]. 

By virtue of the laws of hydraulics, the turbine power 
Pt is calculated by the formula: 


Prit) = h{t)Q{t), 

where h is a pressure drop [Pa = ^^], Q is the water 
flow through the turbine [^]. 

The water flow through the turbine is related with 
h and the opening (closing) guide vanes in the following 
way: 

Q{t) = 

fi{t) = fio + 

where /i is a position of guide vanes [p.u.], /io is a 
given position of guide vanes [p.u.], jiA is the devia¬ 
tion of given position of guide vanes [p.u.], C = S/s/p 
is a constant depending on the construction of penstock 
[ ™v/^ ], S' is a sectional area of water conduit [m^], p 

is the density of water [^]. 

The differential equation, describing the change of 
water flow through the turbine, can be written accord¬ 
ing to the Newton equation for the liquid column en¬ 
closed in the penstock. Then, taking into account that 
S is a sectional area of penstock [m^], I is a length of 
penstock [m], p is a density of water [^], Pu is a con¬ 
stant pressure on the upper end of penstock [Pa], pi is 
a constant pressure on the lower end of penstock (after 
turbine) [Pa], one obtains 

S 

ip 
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Thus, the dynamics of hydraulic turbine is described 
by the system of differential equations 


J UJ = 


Q- l^[Pu 


k 




C^iHo + w 


■Pi 


-Mg, 


C^ipo + paY 


where Mq is a generator torque in the physical values. 

In order to pass to the rotating coordinates (c^, g), 
it is necessary to make the following transformation 



Fig. 6 Scheme of automatic speed governor of hydraulic tur¬ 
bine 


UJA{t) = Uj{t) - CJQ. 


( 5 ) 


Since for the description of dynamics of synchronous 
generator there are used equations in per unit values, 
we write the motion equation of turbine also in per unit 
values: 




Tj \C^{po + piAYu:l{l + s) 


Mg^ , 


where _ 

k 

and in Section there was dehned the generator torque 
in per unit values 

Mg = - ^qid- 

Note that the obtained system contains one control 
signal which corresponds to automatic speed gov¬ 
ernor of turbine. 


5 Mathematical model of speed governor 

Let us consider the simplihed scheme of automatic speed 
governor of turbine [33l pp. 153-157], presented in Fig.[^ 
In this work parameters of hydropower unit are chosen 
so that the speed governor does not get to saturation. 
To simulate the case of saturation it is necessary to 
consider more complex models of speed governors (see, 
e.g., [ 6 l[ 3 l[ 37 l| 3 a| 39 ]). 

The governor has the deadband z, which is speci- 
hed by the technical conditions. The deadband of the 
hydraulic turbine is 30 mHz [33] . The deadband in the 
per unit values is z = 0, 002 . 

The input signal is a relative deviation of rated an¬ 
gular speed s. After the input signal passes through the 
deadband, one obtains the signal 77 ^, which corresponds 
to a signal of measuring device: 


where a is a transmission coefficient of open-cycle con¬ 
trol system. 

The control signal 77 is formed by the formula 


V = -Vs - 


where JIa is a signal of rigid negative feedback. 

Then the control signal is cut off because of restric¬ 
tion on the velocity of change of guide vanes position: 


{ ^, Po<V < Pc^ 

Po, V < Po, 

Pc, Pc < V, 

where po, pc are the maximum velocities of opening and 
closing the vanes. 

A servomotor is presented by the integrator with 
time constant Tc [s]. The output value of servomotor is 
the relative displacement of guide vanes Jia • The stroke 
of servomotor has limit stops Pmin^ Pmaxi which cor¬ 
respond to the minimum and maximum power of the 
turbine [33]. If guide vanes reach limit stops, at which 
further displacement of vanes in the same direction is 
not possible, then displacement of vanes is stopped. In 
some research works it is recommended to model such 
saturation through direct feedback p [38ll39[l43] . 

Consequently, the control signal takes the following 
form 

C = P- ^p{PA,Po), 

where 


ip{PA,Po) 


P, PA+ Po< Pmin and /9 < 0, 
P, PA + Po> Pmax and p>0, 
0 , otherwise. 


The equation of servomotor motion is as follows 

~ _P-ip{PA,Po) 

Pa — rp : 

c 

where 

Pa = Xp.{PA,Po) = 


( a{s — zl2\ s > ^/ 2 , 
Va = <^Xs{s) = < ( t(s + z/2), s < -zjl, 
[ 0 , \s\<zl2. 


pA, pmin PO — PA — Pmax PO, 

— ^ pmin PO^ Po T PA ^ pmin-) 

Pmax P0:i PO pA P‘ 
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In other words, the signal {ia is the signal passed 
through the saturation. Consequently, equation §) can 
be rewritten in the form 

_ P ~ 

T, ’ 

pmin PO — Pa ^ pmax PO' 

Thus, the automatic speed governor of the turbine 
can be described by the following differential equation 

Xp (-CTXs(s) - XpifJ-A, Mo) - ^p{^^A, Mo)) 

M/i = --• 


6 Mathematical model of hydropower unit. 
Operating modes 

Using equations of each structural element of hydropower 
unit, one writes the equations of hydropower unit with 
automatic speed governor in per unit values: 

Oa = ^0 <§, 

T;(c2(/.o + MA)2‘^o'a + s) 

^ CH^io + fiAv)' 

Wd = -a;o(l + s)^q -uiQrid + uioU sin(0o + Oa), 

<Pq = wo(l + s)<l/d - ojQriq - uqU cos( 6 'o + Oa), 

Tr^r = Er- Eg, 

^rd ^rd Erq, 

Erg ^rq ~ Ej-d, (7) 

Xp (-o-Xs(s) - Mo) - ^p{^^A, Mo)) 

M/i = --> 

~ '^d T Eq T ^rq-) 

^q — ^q '^q ^rdi 

xir! . ^ Xod ^ 

Eq -{ - Erq^ 

Xr Xr 

2 

^rd = ^id + ^.9 + — {Eg + Erg) , 

Xrd Xrd 

X^ T 

i7/ — 

^rq — ^q X^rd X^rd- 

Xfq Xrpq 

In the case of the steady-state stability of hydropower 
unit the turbine rotates with a constant angular speed 
{s = 0) and the generator produces the constant power. 
Such an operation mode of hydropower unit corresponds 
to a steady-state (operating) mode of power network. 
Note that for computations the hydropower unit are 
often presented as a source of current or electromotive 
force, the power of which corresponds to a power of 
generator. 


Dynamical stability of hydropower units is consid¬ 
ered in terms of maintaining a given mode. This means 
that if sudden, significant changes of network mode 
arises, then after the transient processes the output 
power of hydropower unit must correspond to the re¬ 
quired power. For example, after short circuits in one 
or more power lines, blackouts, changes of the external 
load, etc. Note that the hydropower unit in the con¬ 
sidered processes can not be represented as a source 
of current or electromotive force since dynamical pro¬ 
cesses in this case have a significant effect both on a 
hydropower unit, and on the network mode. 

An operating mode of hydropower unit corresponds 
to the asymptotically stable equilibrium point of system 
0 . The equilibrium points are the following 


s®* = 0 (i.e. = 0 ), m 2 = 0 , 

6»2 = const, = CmoVPu - Pi, 


TT'St 77' ZTiSt n fpSt (A 

L/g — -t^rq — ^5 -t^rd — ^5 

J-(_21[/sin6'-f/cos6» + i;^), 

H + XdXg Xq 

c = —- {- — UsinO -UcosO- EP, 

^ -h XdXq r 


^d — Xd^d T 

Xrd Xj^d 


rp j 

q -Eqiq^ 


^st _ 


—id 


• Ej>^ 


dyst _ ^ad ■ I ^ 

^rd — ~—^d + y— Er, 


K = 

^ Xrd 


la. 


The position of guide vanes /io is defined from the fol¬ 
lowing equation, which will be called the balance equa¬ 
tion between the turbine and the generator torques: 


Mt{po,(Ao) — Mg{U,0o), 


( 8 ) 


where 

AX f ^ ^E{pu-Pl)^ 

Mt[Po,^o) — - 2 - 

LUq 

Mg{U, 0) = Eil—^Gu^Xd Sin^ 0 + U^Xg cos^ 0- 

ExdXqY 

r‘^-XdXq 2 . r.XdXq-r‘^^^ . 

- -U sm 6 'cos 6 ' — Er - - - U sm^+ 

^ ^ r 

E2xqErU cosO) -z- - - i—XdU sinO -\-rU cos 6 >)+ 

-h XdXq 

^ rXqErjXd- Xq) _ vE^ 

(r^ -h XdXqY (r^ -h XqXd) ' 

Graphical solutions of the balance equation depend¬ 
ing on voltage change in the power network are pre¬ 
sented in Fig. The parameter corresponding to the 
voltage is represented as 1 / = yUnom, where 7 > 0 . 

Equation 0 contains the input parameter U and 
the variable /io- The equality of turbine and genera¬ 
tor torques is attained due to the variable /io- Recall 
that Clio = 14,954rad/s, ^0 = arccos(0.9). Then /io, de¬ 
pending on the voltage U, is found from the balance 
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Fig. 7 Graphical solution of the balance equation depending 
on U. Lines 1, 1’ for 7 = 1 , lines 2, 2’ for 7 = 0.9, lines 3, 3’ 
for 7 = 0.8, lines 4, 4’ for 7 = 0.1 


equation ^ of torques of turbine and generator. A 
plot of /io against U is presented in Fig. 



Fig. 8 Plot of /To against U, k = 40, Oq = arccos(0.9) 

For sufficiently large and small values of [/, the pa¬ 
rameter /io falls on the saturation that corresponds to 
the limit position of guide vanes. In particular this mode 
corresponds to the start of hydropower unit. Then it is 
necessary to consider a more complex balance equation 

MT{fJ^o{U),LUo) = Mg{U,0o + Oa), 

and the balance will be achieved with help of Oaj i-e. 0^2 
may not be equal to zero. In Fig.j^a plot of ^ = ^0 + 
against voltage U is shown. 

If /io does not fall on the saturation, then = 2'Kk. 
Otherwise, = C + 27r/c, where C = const. 

Since for /io = jJ^max (/^o = Mmin) there are restric¬ 
tions on one can get again the saturation. In this case 
the balance will be achieved with the help of uja- 

However for the considered allowed voltage U the 
value jiio{U) does not fall on the saturation, thus we do 
not consider the cases of achieving the balance with the 
help of Oa and 00 a- 



Fig. 9 Plot of 0 against U 


The instantaneous power is determined by the for¬ 
mula 

P{U,e{t)) = {id{t)U sm{d{t)) + ig{t)U cos{d{t))). 

Represent the instantaneous power in the following form 

p{u,e{t)) = Po + PA{u,d{t)), 

where Pq is the required (nominal) power, which is de¬ 
termined by the formula Pq = P{U, 6>o), Pa is a devia¬ 
tion of the required power. 

The change of U leads to the change of power. A 
plot of P against U is shown in Fig. 



Fig. 10 A plot of P{U,0o) 


7 Calculation of generator and turbine 
parameters for the Sayano-Shushenskaya 
hydropower plant 

The radial-axial vertical hydraulic turbines RO-230/833- 
B-677, connected with synchronous generator on the 
umbrella type SVF-1285/275-42 UHL4, are installed at 
the Sayano-Shushenskaya hydropower plant. 

The rest of the system parameters are determined 
as follows: 
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Table 1 Parameters of synchronous generators (SVF- 
1285/275-42 UHL4) [13] 


Parameter 

Value 

Rated angular 
speed, ujQ 

142.8 [rad/s] 

Stator resistance, r 

0.0034 [p.u.] 

Leakage inductive reactance 
of stator winding, Xg 

0.184 [[p.u.] 

Synchronous inductance 
along the axis d, Xd 

1.58 [p.u.] 

Synchronous inductance 
along the axis g, Xq 

0.97 [p.u.] 

Transient resistance, 

0.43 [p.u.] 

Sub-transient reactance 
along the axis d, x'^ 

0.3 [p.u.] 

Sub-transient reactance 
along the axis q, x'^' 

0.31 [p.u.] 

Time constant 
of field winding, Tr 

8.21 [s] 

Sub-transient time 
of field winding, T'J 

0.143 [s] 

Moment of inertia, J 

25.5 • 10^ [kg • m^] 

Field voltage, Er 

530 [V] 

Stator voltage, U 

15.75 [kV] 


Xd = 1.58 [p.u.], Xq = 0.97 [p.u.], Tr = 8.21 [s], J = 
25.5-10® [kg-m^], Er = 530 [p.u.], C = 0.27 [m^^/ykg], 
Xad = 1.396 [p.u.], Xaq = 0.786 [p.u.], Xr = 1.6946 [p.u.], 
Xrd = 1.6155 [p.u.], Xrq = 0.9361 [p.u.], Trd = 0.8666 
[s], Trq = 0.7604 [s], S = 7r/4 ■ 7.5^ [m^], I = 192 [m], 
p = 0.98 ■ 10®[kg/m^], pu = 2.7 ■ 10® [Pa], pi = 0.35 ■ 10® 
[Pa], fc = 40 [kg ■ rv?/s% Qniax = 358 [m®/s]. 

8 Local analysis 

Let us study the local stability of equilibrium points 
of system 0 - It is enough to carry out the analysis 
of stability on the interval [0, 27r) since the solutions 
of the system are 27r-periodic. The equilibrium points 
with respect to 0 are defined from balance equation 0- 
System 0 may have 0, 1, 2, 3 or 4 equilibrium points 
on the interval [0, 27r). 

Let us find the Jacobian matrix of the right-hand 
side of system 0 . For this purpose algebraic system of 
equations ([^ is solved for id^ iq^ Eq^ Erd^ Erq. 


— Impedances of stator winding along the axis d: Xad = 
Xd — Xg = 1.396, 

— Impedances of stator winding along the axis q: Xaq = 
Xa — Xs = 0.786, 


— Impedances of field winding: Xr 


1.6946, 


Xd-x'^ 

— Impedances of damper winding along the axis d\ 

-1 


^rd — ^ad 




-^- ) = 1.6155, 

Xad X I 


a) 


— Impedances of damper winding along the axis q: 


+ 


_^ 

\x'^-Xs Xaq J 


0.9361, 

-Xad = 0.2986, 


'^rq '^aq 

— Resistance of field winding: Xgr = 

— Resistance of damper winding along the axis d: Vrd = 

iXrdXd—X^^)Xrd _ Q I 

u^oXdx'^T- -U.i24D, 

— Resistance of damper winding along the axis q: r^q = 

= 0.0823, 

^OXqT^' ’ 

— Time constant of damper winding along the axis d: 

= 0 . 8666 , 

— Time constant of damper winding along the axis q: 


T — 


^0 P rq 


= 0.7604. 


Table 2 Parameters of turbine (RO-230/833-B-677) 


Parameter 

Value 

Length of penstock, 1 

212 [m] 

Cross sectional diameter of penstock, D 

7.5 [m] 

Available water fiow through turbine, Qnom 

358 [m^/s] 


Thus, for modeling the parameters of hydropower 
unit of the Sayano-Shushenskaya hydropower plant were 
used (HSj): ujq = 27rI42.8/60 [rad/s], r = 0.0034 [p.u.]. 


id = XdEd - XrEr + Xrd^rd. iq = YqEq - Yrq^rq, 


Eq — ZdEd Zrpd^Erpd^ Ej-d — Eq^q ErqErqi 

Erq — QdEd Y Qr^r Y QrdErd: 

where 




a4 ae — I 


h 


Yq = 


_a6-I _I-a4 

yvr — , 1 y^rd — 


bi ’ 


bi ’ 


_ - _ Y — _ - _ 

^ 1 rq — , 

a2 — (27 (22 — (27 


^ _ (6i - (ai - as) (a4 - I)) ^ 

5i(I-a6) ’ ’ 

hi Y (ai — as) (I — a4) 


Erd — 


Eq = 


hi (I — as) 

(27 ^ Ci2 


p — 

1 rq — 


Qd 


<22 — dl d2 — Cij 

(«i - as) («4 ae - 1) - bi 

{bi (1 - a4) 


_ (fli - as) {ae - 1) - 6i _ as - ai 

L/r L/i \ 1 ^rd 


bi (1 - a4) 


''ad 


QjI — Xd-) (22 — Xq^ CL^ — , a4 

Xr Xr 


hi ’ 
Xad 


as 


±gA 

Xrd 


ae 


Xad 

Xrd 


aj — 


^aq 

5 

^rq 


bi = (ui - ae) (a4 - 1) - (ui - as) (1 - ae). 
Then nonzero elements of the Jacobi matrix 

J = {ii,/c}/c=1...9 
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of the right-hand side of system 0 in stationary point 
are defined by the formulas: 


jl,2 = ^0, 


j2,2 = 


TJC^uJl^^r 


i2,3 = 




J2,4 = 


{Xd -Yg)^g+ Yrg 'Prq 


J2,5 = 


-Yg 'I'd + Xd'I'd- Xr 'I'r + Xrd %d 


-Xrl'q . Xrd^q . YrqI'd 

J2,6 — Tf; -; J2,7 — — Tf, -; J2,8 — -; 

J-J J-J J-J 

. _ 2fc(g**)3 . _ 25g** 

2 S' 

is,9 = - , ^ i , j4,i=u;oUcos(ffo+ff^J), 

i pu Pq 

j4,2 = -Wo 'Pq, j4,4 = “Wq r Xd, j4,5 = “Wq 


2 

1.5 


0 


-f++++++++++++++++++++++++++++++++++++++++++++i 


0.4 0.5 0.6 0.7 0.8 


>«+++++++++++++++++++++++++++++++++++++++ 

Aq i Li 1:2 y 


Fig. 11 Stability of equilibrium states of system ([^, defined 
by the least-squares method: blue pluses are stable equilib¬ 
rium states, red crosses are unstable equilibrium states 


9 Analysis of transient processes 


J4,6 — ^0 ^ J4,7 — ^0 ^ ^rdi 

75,1 = Wo C sin( 6>0 + 6*2i)> 75,2 = wq Pd, 

75,4 = Wo, 75,5 = —Wo r Yg, = wo r Yrg, 

Zd . . Z^d 

76,4 = - 77^5 76,6 = -T^T) 76,7 = 

_L rp J. p J. p 

Q d . Qr . Qrd 

J7,4 — 7^, — j7,7 — -7f^^ 

J-rd J-rd J-rd 

P P 1 1 

J8,5 = 7^, J8,8 = J9,2 = 7^, J9,9 = 

Irq Irq 1 c 1 c 

In order to study the local stability of the equilib¬ 
rium states of system 0 relative to the voltage, the 
standard function Isqnonlin of the application package 
MatLab was usecQThis function is based on the least- 
squares method, i.e. on an iterative approximation to 
equilibrium state, that allows one to reduce the com¬ 
putational error and to define more exactly the inter¬ 
val of instability of the system. The initial data for 
this method is putative equilibrium state. Using this 
method the Jacobi matrix calculated at an equilibrium 
is found, and then the Routh-Hurwitz stability crite¬ 
rion is applied. Results are presented in Fig. El the 
equilibrium state 0 = Oq \8 unstable for 7 G [71 ~ 
0.86; 72 ~ 0.9] and stable for other 7, the rest three 
equilibria states are always unstable. 

^ First, the study of local stability of equilibria was carried 
out via eigeuvalues of the Jacobiau matrix aud the staudard 
fuuctiou eig from MatLab. It turued out that all equilibrium 
states are uustable for all 7. Further, the stability criteriou 
of Routh-Hurwitz for the Jacobiau matrix was applied. It 
showed that the equilibrium state 0 is uustable for 7 G 

[71 ~ 0 . 84 ; 72 ~ 1 . 17 ] aud stable for other 7. The obtaiued 
results of these methods cau be explaiued by the fact that all 
calculatious are made with some uumerical error aud because 
of the substautial differeuce iu maguitude. 


During the operation of hydropower unit the transient 
processes related to sudden changes of the work pa¬ 
rameters of the hydraulic unit often occur. As a re¬ 
sult, the following problem arises: to find parameters, 
under which the hydropower unit pulls in the new op¬ 
erating mode after transient processes. This problem 
is closely related to the limit (ultimate) load problem, 
which arises in practice of operation of electrical mo¬ 
tors after sudden change of load torque on the shaft [3l 
[4^[T^[8l [4l[ ^ . For its solution the equal-area method 
is widely used in engineering practice. This method 
was used for some models in the works of A. A. Yanko- 
Trinitskii [42]. In our work modern methods of numer¬ 
ical integration of the system (Runge-Kutta method) 
are combined with the analysis in the spirit of the clas¬ 
sical ideas of Yanko-Trinitskii. 

First the numerical analysis of transient processes 
was carried out with the initial data taken from a small 
neighborhood of the equilibrium state. It is verified 
whether the trajectory goes out from this neighborhood 
after a long integration time (1000 s) or not. As a re¬ 
sult of the study it was obtained that the equilibrium 
state 0 = corresponding to the operating mode, is 
unstable for 7 G [71 ~ 0.85; 72 ~ 0.91]. The received in¬ 
terval of instability is consistent with the local analysis 
and corresponds to the interval [Si^ S 2 ] in Fig. It 
can be defined more exactly due to coefficient k, corre¬ 
sponding the turbine used at the Sayano-Shushenskaya 
hydropower plant. 

Further, the numerical analysis of transient processes 
was carried out with various initial data. For values 7, 
corresponding the local stability, hidden oscillation^ 

^ An oscillation can generally be easily numerically local¬ 
ized if the initial data from its open neighborhood in the phase 
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are not found numerically and all simulated trajecto¬ 
ries attract to the equilibrium states. For values 7, cor¬ 
responding to the local instability, simulated trajecto¬ 
ries attract to self-excited periodic solutions. The local 
bifurcation, in which an equilibrium loses stability and 
a small stable limit cycle is born, occurs in considered 
multidimensional system (this bifurcation is an analog 
of Andronov-Hopf bifurcation miniEi]). 

According to [36], immediately before the accident 
the power of the second hydropower unit was 475 MW 
at a head of 212 meters (i.e., it worked in the not rec¬ 
ommended zone II (Fig.[^. On the day of the accident 
the power of the second hydropower unit was reduced in 
accordance with the commands of the group controller 
of active and reactive power. 

Below three cases are modeled: 

1. operation of hydropower unit at the rated voltage 
(Fig. 0 point A) with initial data 

( 0 , 1 , 0 , 0 , 0 , 0 , 0 , 0 , 0 ), 

2. reducing the power of hydropower unit, that cor¬ 
responds to reducing voltage to 0.89 of the rated 
voltage (Fig. point B), 

3. reducing power of hydropower unit, that corresponds 
to reducing voltage to 0.7 of the rated voltage (Fig.[^ 
point C). 

The results of modeling have shown that at the 
rated voltage the trajectory of the system after tran¬ 
sient processes is attracted to the equilibrium state, 
which corresponds to operating mode of hydropower 
unit (Figs. [T^ . Further at some instant the volt¬ 
age is reduced to 0.89 of the rated voltage. In this case 


space lead to a long-term behavior that approaches the os¬ 
cillation. Therefore, from a computational perspective, it is 
natural to suggest the following classification of attractors 
|2III3III32II28III9| . which is based on the simplicity of finding 
their basins of attraction in the phase space: An attractor is 
called a self-excited attractor if its basin of attraction inter¬ 
sects with any open neighborhood of an equilibrium, otherwise 
it is called a hidden attractor. For a self-excited attraetor 
its basin of attraction is connected with an unstable equi¬ 
librium and, therefore, (standard eomputational proeedure) 
self-excited attractors can be localized numerically by the 
standard eomputational proeedure: by constructing a solution 
using initial data from an unstable manifold in a neighbor¬ 
hood of an unstable equilibrium, and observing how it is at¬ 
tracted, and visualizing the oscillation. In contrast, the basin 
of attraction for a hidden attractor is not connected with any 
equilibrium. For example, hidden attractors are attractors in 
systems with no equilibria or with only one stable equilibrium 
(a special case of the multistability: coexistence of attractors 
in multistable systems). Well known examples of the hidden 
oscillations are nested limit cycles in I6th Hilbert problem 
(see, e.g., |20II28| ) and counterexamples to the Aizerman and 
Kalman conjectures on the absolute stability of nonlinear con¬ 
trol systems |24ll27ll7ll28| . Hidden oscillations in the models 
of electrical machines are discussed, e.g., in 


the trajectory of the system after transient processes 
is attracted to the stable limit cycle (Figs. [T^, he., 
vibrations are arisen in the hydropower units. Then the 
voltage is reduced to 0.7 of the rated voltage. The tra¬ 
jectory of the system after transient processes is at¬ 
tracted to equilibrium state (Figs.j^, 17). 


, , , , , , 

0^ -deviation of operating angle 


s -relative frequency deviation 



- deviation of given vanes position 



-deviation of given power 


Fig. 12 Stable equilibrium in the mathematical model of 
hydropower unit, U = 15.75 • 10^ [V] 


stable equilibrium state 



Ga -deviation of 
operating angle 


s - relative frequency deviation 

Fig. 13 Stable equilibrium in the mathematical model of 
hydropower unit - projection onto (0z\, s), U = I5.75T0^ 

[V] 

































A dynamical model of Sayano-Shushenskaya hydropower plant: stability, oscillations, and accident 


11 



-deviation of given power 


Fig. 14 Limit cycle in the mathematical model of hy¬ 
dropower unit, U = 0.89 • 15.75 ■ 10^ [V] 



Fig. 17 Stable equilibrium in the mathematical model of 
hydropower unit - projection onto s), U = 0.7-15.75- 

10^ [V] 


From Fig. it is clear that the maximum value of 
the oscillation amplitude is reached at 7/ = 0.897/nom- 



Fig. 15 Limit cycle in the mathematical model of hy¬ 
dropower unit - projection onto (0^, s), U = 0.89-15.75- 

10^ [V] 



Fig. 18 Amplitude of oscillations for different voltages U = 
PUnom (/5 - percent of capacity) 

Thus, the results of modeling are sufficiently consis¬ 
tent with the full-scale tests carried out for hydropower 
units of the Sayano-Shushenskaya hydropower plant. 


^..■■■■■ ■■ ■ 

0^ -deviation of operating angle 

—- - --- 

”550 555 560 565 570 575 580 585 590 

s -relative frequency deviation 

o^A/WWVWVV\/|^ . .I ^ ■ 

\x^ -deviation of given vanes position 

. — — - . I — 

-deviation of given power 


595 600 


Fig. 16 Stable equilibrium in the mathematical model of 
hydropower unit, U = 0.7 - 15.75 - 10^ [V] 
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